Quantifying non-Gaussianity for quantum information 



Marco G. Genon{3 
CNISM, UdR Milano, 1-20133 Milano, Italia 

Matteo G. A. ParifQ 
Dipartimento di Fisica, Umversita degli Studi di Milano, 1-20133 Milano, Italia 
(Dated: December 13, 2010) 

We address the quantification of non-Gaussianity of states and operations in continuous-variable 
systems and its use in quantum information. We start by illustrating in details the properties and the 
relationships of two recently proposed measures of non-Gaussianity based on the Hilbert-Schmidt 
(HS) distance and the quantum relative entropy (QRE) between the state under examination and a 
reference Gaussian state. We then evaluate the non-Gaussianities of several families of non-Gaussian 
quantum states and show that the two measures have the same basic properties and also share the 
same qualitative behaviour on most of the examples taken into account. However, we also show 
that they introduce a different relation of order, i.e. they are not strictly monotone each other. We 
exploit the non-Gaussianity measures for states in order to introduce a measure of non-Gaussianity 
for quantum operations, to assess Gaussification and de-Gaussification protocols, and to investi- 
gate in details the role played by non-Gaussianity in entanglement distillation protocols. Besides, 
we exploit the QRE-based non-Gaussianity measure to provide new insight on the extremality of 
Gaussian states for some entropic quantities such as conditional entropy, mutual information and 
the Holevo bound. We also deal with parameter estimation and present a theorem connecting the 
QRE non-Gaussianity to the quantum Fisher information. Finally, since evaluation of the QRE 
non-Gaussianity measure requires the knowledge of the full density matrix, we derive some experi- 
mentally friendly lower bounds to non-Gaussianity for some class of states and by considering the 
possibility to perform on the states only certain efficient or inefficient measurements. 
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I. INTRODUCTION 

In the recent years we have witnessed a big effort in the 
theoretical and experimental investigation of continuous- 
variable (CV) quantum information. Gaussian states are 
experimentally produced with an high degree of con- 
trol, especially in quantum optics, and Gaussian mea- 
surements may be effectively implemented in different 
settings. Besides, despite they belong to an infinite- 
dimensional Hilbert space, Gaussian states are easy to 
handle from the theoretical point of view, being fully de- 
scribed by the first and second moments of the canonical 
operators [H-Q- The remarkable role of Gaussian states 
has been highlighted in [4| , where it has been proved that 
they are extremal at fixed covariance matrix for several 
relevant quantities as channel capacities and entangle- 
ment measures and also in the framework of CV quantum 
key distribution in 0-0] , where it has been shown that 
Gaussian attacks are optimal against all individual and 
collective eavesdropping strategies. For these reasons, 
Gaussian states played a prominent role in the devel- 
opment of CV quantum information and, as a matter of 
fact, most of the protocols designed for finite-dimensional 
Hilbert spaces have been firstly translated in the CV set- 
ting for Gaussian states p|. 
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In the recent years, however, it has been realized that 
there are situations wherein non-Gaussianity (nG) in the 
form of non-Gaussian states or non-Gaussian operations 
is either required or desirable to achieve some relevant 
tasks in quantum information processing. As for ex- 
ample, it is known that nG is crucial for the realiza- 
tion of ent ang lement distillation [UJnJ, quantum error 
correction [12J and cluster states quantum computation 
(l3l ITij ]. Besides, a non-Gaussian measurement and/or 
non-Gaussian states arc crucial to observe violation of 
loophole free Bell tests with continuous variables [15l424| . 
In addition, improvement of quantum teleportation and 
quantum cloning of coherent states can be obtained by 
using respectively non-Gaussian states or non-Gaussian 
operations [25l428| . In turn, bipartite Gaussian states 
have minimum entanglement for given second moments 
and this influences their performances in quantum in- 
formation protocols. Non-Gaussian operations also find 
application in noiseless amplification [29l . |30j obtained 
in conditional fashion, whereas non-Gaussian states have 
been proven useful to improve parameter estimation in 
quantum optics [3l|, [32| . The current state of the art is 
schematically depicted in Fig. [1] 

For the reasons outlined above several protocols have 
been designed theoretically [33l . l34j and experimentally 
realized [35-45] to produce single mode or two- mode non- 
Gaussian states, in different physical settings, and in par- 
ticular to perform squeezin g pu rification [461 ] and CV en- 
tanglement distillation (40l T l47l449| . Basically, they may 
be divided into two main categories: those based on non- 
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linear interaction of order higher than two [56J, l57j , as for 
example the Kerr effect [58 60]), and those based on con- 
ditional measurements. Indeed, the nonlinear dynamics 
induced by conditional measurements has been analyzed 
for a large variety of schemes [6lT[75l |. also including, be- 
sides photon addition and subtraction schemes, optical 
state truncation of coherent states fijll. state filtering 
by active cavities [65l . |66| , synthesis of arbitrary unitary 
operators [67} and generation of optical qubit by condi- 
tional interferometry (68[. Conditional state gen eration 
has been achieved in the low energy regime |45ll49l.l50| by 
using si ngle-ph oton detectors, and also in the mesoscopic 
domain [5ll453| . Realisations of non-Gaussian states have 
been also reported in optical cavities [Hj], and in super- 
conducting circuits (55|. 

Being recognized as a resource for CV quantum infor- 
mation the need of quantifying the nG character of states 
and operations naturally arises and different measures of 
non-Gaussianity have been proposed [76T - I781 ]. These mea- 
sures have been used to assess the role of nG in different 
quantum information and communication tasks as tele- 
portation [79j , quantum estimation [31| , experimental en- 
tanglement quantification (80j and entanglement transfer 
between CV states and qubits [U |82[ . In [83| the rela- 
tionship between nG and the Hudson's theorem [84| have 
been studied, obtaining at fixed purity an upper bound 
for non-Gaussian states having a positive Wigner func- 
tion, while in [85| nG bounded uncertainty relations are 
derived. The entropic measure proposed in [77| has been 
used to quantify exactly the nG of experimentally pro- 
duced photon-added coherent states and a lower bound 
has been evaluated experimentally in 86] for conditional 
states obtained via an inefficient photo-detection on clas- 
sically correlated thermal beams. 

In this paper, we address in details the quantification of 
non-Gaussianity of states and operations in continuous- 
variable systems and analyze its use in quantum informa- 
tion. At first, we review the two measures proposed in 
[76| and by studying in more details their properties 
and their relationships and then exploit them to assess 
some relevant Gaussification and de-Gaussification pro- 
tocols, and to investigate the role of non-Gaussianity in 
entanglement distillation and quantum communication. 

The paper is structured as follows: in the next section 
we briefly review some topics on the quantification of 
non-Gaussianity of a classical probability (density) dis- 
tribution. In Sec. IIIII we introduce notation and give the 
definition of Gaussian states along with their properties, 
while in Sec. IIVI we review the two quantum measures 
of nG, proving their properties and highlighting the re- 
lationships between them. In Sec. |V]we evaluate the nG 
measures for some relevant non- Gaussian states, com- 
paring them and observing if and when they give the 
same order relation. In Sec. I VII we employ the two mea- 
sures of non-Gaussianity to address Gaussification pro- 
cess due to the interaction of the system with a bath of 
harmonic oscillators in the vacuum state (i.e. dissipation 
a zero temperature) and de-Gaussification process due 
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FIG. 1. The departure from the Gaussian world is required 
to achieve specific task such as the distillation of Gaussian 
entanglement, universal quantum computation with Gaussian 
cluster states and the violation of loophole-free Bell tests with 
continuous variables (the lower box) and lead to an increase of 
entanglement at fixed covariance matrix, with improvement 
of relevant protocols such as teleportation and cloning of co- 
herent states and parameter estimation of both unitary and 
lossy channels (the upper box). 



either to phase-diffusion or to self-Kerr interaction. In 
Sec. I VIII we study the role of nG in two paradigmatic 
examples of entanglement distillation protocols, while in 
Sec. I Villi we show how the amount of nG is related 
to some entropic-informational quantities as the Holevo 
bound, conditional entropy and mutual information. In 
Sec. IIXI we deal with parameter estimation and present 
a theorem relating the non-Gaussianity and the quan- 
tum Fisher information. Finally, in Sec. |X] we address 
the experimental evaluation of nG in situations where 
state tomography is not available, and present some ex- 
perimentally friendly bounds for the estimation of nG of 
some classes of states. Sec. IXII closes the paper with 
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some concluding remarks. 



II. NON-GAUSSIANITY OF A CLASSICAL 
PROBABILITY DISTRIBUTION 

According to the central limit theorem the Gaussian 
distribution is ubiquitous in the description of natural 
phenomena. In turn, deviations from the Gaussian be- 
haviour are often the sign that some interesting phe- 
nomenon occurs [90l - [95| . and thus a considerable atten- 
tion has been devoted to the detection and quantification 
of non-Gaussianity of a classical distribution. Basically, 
there are two main approaches: The first one is based on 
the evaluation of higher moments of the distribution, in 
particular the third and the fourth central moments, to 
assess Skeweness and Kurtosis of the distribution in com- 
parison to those of a Gaussian one. The second approach 
is based on the evaluation of the Shannon entropy of the 
distribution, upon the fact that Gaussian distributions 
maxmize it at fixed variance. More recently, it turned 
out that non-Gaussianity is relevant in the framework of 
Independent Component Analysis (ICA) [96]. ICA is a 
method developed in the last decades in which the goal 
is to find a linear representation of non-Gaussian data 
so that the components are statistically independent. In 
this method the solution is obtained by the maximisation 
of the nG of the components and thus, to accomplish this 
goal, different measures of nG have been proposed. 

Let us consider a scalar- valued random variable Y with 
a probability density function p(y) := P(Y — y). Its fc-th 
central moments are defined as 



where 



E[{Y-ii) k 
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(1) 



(2) 



is the mean value of the distribution. We say that Y is 
Gaussian distributed if its probability density function is 
a Gaussian function, that is 



p(y) 
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where a 2 = E[(Y — fi) 2 ] is the second moment of the 
distribution, called variance. In the following we will 
present the definition of two measures of nG for a classical 
scalar-valued random variable. 



A. Kurtosis 

The fist ever considered measure of nG has been the 
Kurtosis, that is the fourth-order cumulant defined by 
the formula 



K{Y) = E[(Y - M ) 4 ] - 3a 2 . 



(4) 



K(Y) is zero for a Gaussian random variable while for 
most (but not all) non-Gaussian random variables, takes 
values different from zero. Kurtosis can be positive or 
negative. Random variables that have negative Kurto- 
sis are called sub- Gaussian and are characterized by a 
probability density function with heavy tails. On the 
other hand, the ones with positive Kurtosis are called 
super-Gaussian and they have tipically a "flat" distribu- 
tion (constant near the mean value and very small for 
"distant" values of the variable) . Anyway, typically the 
nG is measured by the absolute value or the square of the 
Kurtosis. While it is relatively simple to evaluate it if the 
probability density function is known, Kurtosis presents 
some drawbacks when its value has to be estimated from 
a measured sample. Indeed, it can be very sensitive to 
outliers and its value may strongly depend on only a few 
observations in the tails of the distribution. For these 
reasons Kurtosis is not considered a robust measure of 
nG. 



B. Negentropy 

Given a discrete random variable X — {x%, . . . , x n } 
with a probability distribution p(xi) = P(X = Xi), we 
can define its (Shannon) entropy as 



H(X) = -Y,P{x i )\ogp{x i ). 



(5) 



This definition can be generalized for a continuous- valued 
random variable Y, in which case it is called differential 
entropy: 



H(Y) = - Jdyp{y) log p(y) . 



(6) 



A fundamental result of information theory states that 
at fixed variance, Gaussian variables have the largest en- 
tropy. Following this result, one may define a measure of 
nG, called negentropy, as 



N{Y) = H(G) - H(Y) 



(7) 



where G is the Gaussian random variable with the same 
variance of Y. Due to the above mentioned result, ne- 
gentropy is always non-negative and it is equal to zero 
only for Gaussian random variables. Negentropy is thus 
well justified by statistical theory but its computation is 
typically very difficult. However, simpler approximations 
based on evaluations of moments of the random variable 



have been introduced and used for ICA purposes [97| • In 
the following we will see that one of the two quantum 
measures that will be analyzed, even if defined starting 
from an another quantity, will result to be the quantum 
analogue of the negentropy here presented. 



4 



III. QUANTUM GAUSSIAN STATES 

For concreteness, we will use here the quantum optical 
terminology of modes carrying photons, but our approach 
may be equally applied to any bosonic (CV) system. Let 
us consider a system of n modes described by mode op- 
erators dfc, k = 1 . . . n, satisfying the commutation rela- 
tions [<Zfe,aj] = Skj- A quantum state g of the n modes 
is fully described by its characteristic function [87} 

X [g](X) = Tx[gD(X)} 

where D(X) = l -Dfc(Afc) is the n-mode displacement 
operator, with A = (Ai, . . . , A n ) T , X k G <D, and where 

D k (\ k ) = exp{Afce4 - \* k a k } 

is the single-mode displacement operator. Analogously, 
quantum states can be fully described by the Wigner 
function, i.e. the Fourier transform of the characteristic 
function: 



W[g](a) 



d 2n X 



A* CK + CK* A 



x[q]W 



(8) 



The canonical operators are given by: 

1 / U 
Qk = -j=(a k + a k ), 

with commutation relations given by [qj,Pk] = 
iSj k . Upon introducing the real vector R = 
(qi,Pi, . . . ,q n ,p n ) T , the commutation relations rewrite 
as 



[R k , Rj 



in 



where £l k j are the elements of the symplectic matrix f2 = 
i®k=i a 2> CT2 being the y-Pauli matrix. The covariance 
matrix (CM) er = cr[g] and the vector of mean values 
X = X [g] of a quantum state g are defined as 



X j - ( R j) 

a k 3 = -^({Rk,Rj}) - (Rj)(Rk) 



(9) 



where {A, B} — AB + BA denotes the anti-commutator, 
and (O) — Tr[oO] is the expectation value of the operator 
O. 

A quantum state go is referred to as a Gaussian state 
if its characteristic function or equivalently the Wigner 
function, have a Gaussian form, in the Cartesian notation 



(10) 



r „ exp \ — ^(Y — X) T cr^ 1 (Y — X)} , s 

W[qq](Y) = 1 2V V '± (11) 

(2 7 r)™ v / Dct[^] V ' 



where A and Y are real vectors, 

A = (ReAi,ImAi, . . . , Re A n , Im \ n f (12) 
Y = (Rc«i,Imai, . . . , Rea„, Ima„) T . (13) 

Of course, once the covariance matrix and the vector of 
mean values are given, a Gaussian state is fully deter- 
mined. As for example, the purity h[qg] = Tr[^c?] °f a 
n-mode Gaussian state may be expressed as 



1 



(14) 



2™Vdet cr 

A n-mode Gaussian state can be always written as 

QG = U s ®l=i v k {n k )ul 
where v k (n k ) = (1 + ri fc ) _1 [n k /(l + n k )] a i ak is a single- 



a\a k v k {n k ) 



average 



mode thermal state with n k = Tr 

number of photons, and Us denotes the unitary evolu- 
tion generated by a generic Hamiltonian at most bilinear 
in the mode operators, that is an evolution correspond- 
ing to a symplectic transformation in the phase-space 
[88j. Any mapping, either unitary or completely-positive, 
transforming Gaussian states into Gaussian states is a 
Gaussian operation. 

For a single-mode system the most general Gaussian 
state can be written as 

g G = D(a)S((Mn)SH()DHa), 

D{a) being the displacement operator and S(Q — 
exp[i£(a'i') 2 — \C,*a 2 ] the single-mode squeezing opera- 



tor with a, C = reU 
matrix has entries 



G (D. The corresponding covariance 
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o~n = (n + [cosh(2r) — sinh(2r) cos(<p)] , (15) 

0-22 = (n+-) [cosh(2r) + sinh(2r) cos(<p)] , (16) 

o~i2 = o~2i =('">+-) sinh(2r) sin(y) . (17) 

The Von-Ncumann entropy S(g) = —Tr[g logg] of a 
single-mode Gaussian states may be written as 



1 



1, 



S(q g ) = h(Vdsta) = h(-) = h(n + -) 
= (n + 1) log(n + 1) - n log(ra) , 
where we have introduced the function 

h{x) = (x + i) log(x + i) - (a; - ^) log(i - ^] 



(18) 



(19) 



For a two-mode Gaussian state, the covariance matrix 
is a real is a real 4x4 symmetric definite positive block 
matrix with ten independent parameters 



(20) 



A 




C T 


B 
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Matrices A, B and C are 2x2 real matrices, represent- 
ing respectively the autocorrelation matrices of the two 
modes and their mutual correlation matrix. Any two- 
mode CM <j may be brought to its standard form local 
symplectic operations, i.e. local Gaussian operations. In 
the standard from, matrices A and B are proportional to 
the identity and C is diagonal. Using the four local sym- 
plectic invariants I\ = det(A), I 2 = det(S), I3 = det(C), 
I4 = det(cr) the symplectic eigenvalues, denoted by d± 
with d_ < d 4. read as follows 



d± = 



where A(cr) = I\ + I 2 + 2/3. Using the symplectic eigen- 
values, the uncertainty relation re- writes as <i_ > 1/2 
and the Von-Neumann entropy as [HI 



S( QG ) = fc(d_) + h(d+) . 



(22) 



IV. 



QUANTUM NG MEASURES: DEFINITIONS 
AND PROPERTIES 



In this section we review the definitions of the nG 
measures for quantum states proposed in [76l |77| and 
illustrate in details their properties. Although the two 
measure are based on different quantities, as the Hilbert- 
Schmidt distance and the quantum relative entropy, they 
share the same basic idea: one wants to quantify the non- 
Gaussianity of a quantum state g in terms of the distin- 
guishability of the state itself from a reference Gaussian 
state r, chosen as the Gaussian state with the same first 
and second moments of g, that is such that 



X[t]=X[q] 
<t[t] = <t[q] . 



(23) 



Notice that a sim ilar line of reasoning has been adopted 
in Refs. f98H100l | to define a measure of non-classicality 
via the Hilbert-Schmidt distance. Here, roughly speak- 
ing, the two nG measures provide the quantization of 
the classical approaches to asses non-Gaussianity based 
on moments and negentropy respectively. In the follow- 
ing we review their properties and also provide a critical 
comparison with another quantities proposed in litera- 
ture [zli- 



A. Measuring the non-Gaussianity using 
Hilbert-Schmidt distance from a Gaussian reference 

Given two quantum states g\ and g 2 , the Hilbert- 
Schmidt distance is defined as 



Vhs[qi,Q2\ = ( ^Tr[(gi - g 2 f] 



1/2 



1/2 



where n[g] is the purity of g and K[g\,g 2 ] = Tr[pi g 2 ] 
denotes the overlap between gi and g 2 . We define the 
degree of non-Gaussianity of the state g as the squared 
renormalized HS distance I7al 



Sa[q] 



iAq\ 



(24) 



of the state g from the state r, which is a reference Gaus- 
sian state chosen as in Eq. (f2"3"|) . The relevant properties 
of <5a[p] are summarized by the following Lemmas: 



(21) Lemma Al 8a[q] = iff g is a Gaussian state. 



Proof: If 6a[b] = then g = r and thus it is a Gaus- 
sian state. If q is a Gaussian state, then it is uniquely 
identified by its first and second moments and thus the 
reference Gaussian state r is given by r = g, which, in 
turn, leads to Dug[g, r] = and thus to <5a[p] = 0. □ 

Lemma A2 If U is a unitary map corresponding to a 
symplectic transformation in the phase space, i. e. if U = 
exp{— iH} with hermitian H and at most bilinear in the 
field operators, then 8a[U g U>] = <5a[£>]- 

Proof: Let us consider g 1 ~ UgU^. Then the covari- 
ance matrix transforms as cr[g'] = S <r[g] E T , S being 
the symplectic transformation associated to U. At the 
same time the vector of mean values simply translates to 
X' = X + Xq, where -X"o is the displacement generated 
by U. Since any Gaussian state is fully characterized by 
its first and second moments, then the reference state 
must necessarily transform as r' = Ut U\ i.e. with 
the same unitary transformation U . Since the Hilbert- 
Schmidt distance and the purity of a quantum state are 
invariant under unitary transformations the lemma is 
proved. □ 

This property ensures that single-mode displacement 
and squeezing operations, as well as two-mode evolutions 
as those induced by a beam splitter or a parametric am- 
plifier, do not change the Gaussian character of a quan- 
tum state. The lemma also allows us to always consider 
state with zero mean values. 

Lemma A3 8a[q] is proportional to the squared L 2 ((D") 
distance between the characteristic functions ( or alter- 
natively the Wigner functions ) of g and of the reference 
Gaussian state r . In formula: 



(25) 
(26) 



8 A \g]<x J d 2n X [ X [g]{\) - X [t](\)Y , 

S a [q]oc Jd 2n a [W[g\{cx)-W[T]{a.)f 

Proof: Using the identities 
d 2n \ 



Tr[Oi0 2 ] 



X [Oi](A)x[0 2 ](-A) 



(27) 



vr" J d 2 ™aIU[0 1 ](a)IU[0 2 ](a) (28) 



6 



and the fact the characteristic functions of self-adjoint 
operators are even functions of A we obtain 

Ks[Q,r] = - / — [xM(A)-xM(A)] 2 (29) 



d' n \{W[e]{oc)-W[r]{oc)Y , (30) 



which proves the Lemma. □. 

Since the notion of Gaussianity of a quantum state 
is connected to the shape of its characteristic (Wigner) 
function, and since the characteristic function of a quan- 
tum state belongs to the L 2 ((D ra ) space [83|, we address 
L 2 ((D) distance to as a good indicator for the non Gaus- 
sian character of g. 



Lemma A4 Consider a bipartite state g = gA® Qg- 
go is a Gaussian state then 5 a [q] — 5a [qa] ■ 

Proof: we have 

lAe] = iAqa]iAqg] 

[i[t] = fi[T A ]fj,[r G ] 
k[q,t] = k[qa,t a ]k[q g , Qg] ■ 

Therefore, since k[qg, go] = iAbg] we arrive at 



If 



Sa[q] = 



iAqa]iAqg] + h[taMqg] - 2n[g A , t A ]k[q g , Qg] 



2n[g A ]n[gG] 



= 5 A [g A ].n 



(31) 



Notice, however, that 5 a [q] is not generally additive (nor 
multiplicative) with respect to the tensor product. If we 
consider a (separable) multi-partite quantum state in the 
product form g = ®JJ =1 £>fc, the non-Gaussianity is given 
by 



Sa[q] 



nLi Hgfc] + nLi mm - 2 nLi K \ek, r k ] 
2 riLiA i bfc] 



( 32 ) 

where Tfe is the Gaussian state with the same moments 
of gk- In fact, since the state g is factorisable, we have 
that the corresponding Gaussian r is a factorisable state 
too. 

For single-mode quantum states we have collected sev- 
eral numerical evidences that 5A[g] = 1/2 represen ts an 
upper bound for the HS nG of any quantum state [lOlj j . 
The same conclusion is indirectly suggested by the results 
obtained in [83[ and this leads to formulate the following 

Conjecture A5 For single-mode quantum states we have 
that 5a [q] < o • 

In particular, the conjecture has been numerically veri- 
fied for single-mode CV states expressed as finite super- 
position of Fock number states, i.e. for truncated states 
of the form g — J2n k=o Qnk\n){k\. We have generated 
at random a large number states for varius values of the 



truncating dimension N and evaluated the correspond- 
ing nG 5a ■ Results have shown that the value of the nG 
5a is bounded by 1/2 and the typical nG (the value of 
5a with the largest occurence) decreases with both the 
purity and the truncating dimension. 



B. Measuring the non-Gaussianity using the 
quantum relative entropy to a reference Gaussian 

Given two quantum states g\ and £>2, the quantum 
relative entropy (QRE) is defined as 



S(£>i||ft2) = Tr[pi(log£>i - log 2 )] 



(33) 



As for its classical counterpart, the Kullback-Leiber di- 
vergence, it can be demonstrated that < S(^i || ^2) < 00 
when it is definite, i.e. when the support of the first state 
in the Hilbert space supp g\ C supp Q2, is contained in 
that of the second one. In particular, S(gi||g2) = iff 
qi = Q2- This quantity, though not defining a proper 
metric in the Hilbert space (it is not simmetric in its 
arguments), has been widely used in different fields of 
quantum information as a meas ure of sta tistical distin- 
guishability for quantum states [l02l Il03j because of its 
nice properties and statistical meaning. In fact if we con- 
sider two quantum states g and r and we suppose to per- 
form iV measurements on g, the probability of confusing 
q with t is (for large N) Pn{q — > r) ~ cxp{— A r S(p||r)}. 

The degree of non-Gaussianity of a state g may be 
quantified as [77| 

5 b [q]=S(q\\t) 

where r is the reference Gaussian state with the same 
first and second moments as in Eq. (j23|) . Notice that, 
because of the choice of r, Tr[rlogr] = Tr[glogr] and 
thus 

5 b [q] = Tr[glogg] - Tr[ e logr] = S(r) - S(g) , (34) 

where S(g) is the von Neumann entropy of a quantum 
state g. The nG measure 5b may be considered as the 
quantum analogue of the negentropy introduced in Sec. 
HTl where the differential entropy is replaced by the von 
Neumann entropy of the quantum states under investi- 
gation. 

At fixed Von-Neumann entropy nG is determined by 
the first two moments of the canonical operators, which 
in turn uniquely determine the reference Gaussian state. 
Upon using formulas from Sect. IHII we may write explicit 
formulas of 5b for single- and two- mode states 



5b[q] — h(Vdetcr) — S(g) single-mode states, (35) 
5b[q] = h(d_) + h(d + ) — S(g) two-mode states, (36) 

where d± are the symplectic eigenvalues of the two-mode 
CM and h(x) is given in Eq. (fi"9"]). 

The relevant properties of 5b\q] are summarized by the 
following lemmas. As a matter of fact, the QRE mea- 
sure of non-Gaussianity owns all the relevant properties 



proved for 8 a, and shows additional properties concern- 
ing the evolution under generic (not unitary) Gaussian 
maps and under tensor product. 

Lemma Bl 8b[q] — iff g is a Gaussian state. 

Proof: If Sb[q] = then g = r and thus it is a Gaus- 
sian state. If q is a Gaussian state, then it is uniquely 
identified by its first and second moments and thus the 
reference Gaussian state t is given by r = g, which, in 
turn, leads to S(g||r) = and thus to <5b[£>] = 0. □ 

Lemma B2 // U is a unitary map corresponding to a 
symplectic transformation in the phase space, i.e. if U = 
exp{— iH} with hermitian H that is at most bilinear in 
the field operators, then 5b[U g W] — Sg[g\. 

This property ensures that single-mode displacement 
and squeezing operations, as well as two-mode evolutions 
as those induced by a beam splitter or a parametric 
amplifier, do not change the Gaussian character of a 
quantum state. The lemma also allows us to always 
consider state with zero mean values. 

Proof: The lemma follows from the invariance of QRE 
under unitary operation. □ 

Lemma B3 5b is additive for factorized states: 5b [gi <8> 
Q2] = <5b[£i] + 5B[g-2\- As a corollary we have that if g-2. 
is a Gaussian state, then 8b[q\ — 

Proof: The overall reference Gaussian state is the tensor 
product of the relative reference Gaussian states of gi and 
Q2 , t = Ti Cg) T2 . The lemma follows from the additivity 
of QRE and the corollary from Lemma BI. □ 

Lemma B4 5b monotonically decreases under partial 
trace, that is, given a bipartite state g, then <5b[pa] < 
5 B [g] and 5 B [gB] < 5 B [g], where g A = Tr B [g], gs = 
Tta[q]- 

Proof Let us consider the partial trace state gA (gs)'- 
its CM is the submatrix of er[g] obtained by dropping 
lines and rows involving expectation values on the sys- 
tem B {A). Analogously, the first moment vector is the 
proper subvector of X[g\. Therefore, the reference Gaus- 
sian state ta (tb) must necessarily satisfies ta = Trs[r], 
where r is the Gaussian reference of g (tb — Tr^r]). 
The QRE monotonically decreases under partial trace 
and thus the lemma is proved. □ 

Actually, the above statement can be strenghtened, as 
it is expressed by the following Lemma. 

Lemma B5 Let us given a generic bipartite state g, we 
have 5[g] > 5[gA] + 5[g B \- 

Proof: It has been shown in jl08| that QRE decreases 
monotonically under a generic (non-linear) coarse grain- 
ing. A simple example of non-linear coarse graining that 
can not be obtained via a completely positive quantum 
map is the operation g — > gA® Qb- Because of this prop- 
erty we have 5[g] > 5[gA <8 £>b] = <5[pa] + 5[gs] where in 
the last equality we have used lemma B3. □ 
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FIG. 2. Non-Gaussianity 5b[iPp\ of PNES states (blue upper 
curves) and of their partial traces 25b [qp] (red lower curves) 
as a function of the total energy of the PNES. Solid curves 
are for TMC, dashed for PSSV and dotted for PASV. 

As an example let us consider the class of pure two- 
mode photon-number entangled states (PNES) 

hM)=£>„|n)|n> (37) 

n 

together with their (equal) partial traces gA and gs , that 
is the diagonal mixtures of Fock states given by 

n 

Relevant examples of non-Gaussian PNES are given by 
the photon subtracted (PSSV) ip n cx (n+l)x n+1 and the 
photon-added two-mode squeezed vacua (PASV) ip n cx 
m"" 1 , which are obtained from the Gaussian PNES 
ip n cx x n < x < 1 (twin-beam state) by the ex- 
perimentally feasible operations of photon subtraction 
p — > a\a2ga\a\ and addition g — > a\a\ga\a2 respectively 
[361 ] ; (iii) the pair-cohe rent or tw o-mode coherently cor- 
related states (TMC) pMpiOrl with Poissonian profile 
ip n cx ^, A € K. The mean energy of PNES is N P = 

MaW + a\aM) = 2N where N = £"= IVgV 
whereas correlations between the modes can be quanti- 
fied by C = ReX^oV'n 1 /Wi( n +l) and entanglement is 
given by the Von-Neumann entropy of the partial traces 
60 = — J2 n V'n 1°§ • m t urn i the covariance matrix of 
a PNES equals that of a symmetric Gaussian state in 
standard form, with diagonal elements equal to + i 
and off-diagonal blocks given by C — diag(C, — C). In 
Fig. [2] we report the nG 5b[iPp] of PNESs as a function 
of the overall energy Np together with the sum of the 
nGs of the partial traces i.e. c$b[£u] + ^s[g6] = 25s[gp}- 
As predicted by the previous lemmas, the nG of a PNES 
state is always larger than the one of its partial traces 
and also of their sum. 
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Lemma B6 8 B [g] monotonically decreases under Gaus- 
sian quantum channels, that is 8[£q(q)] < S[g]. 

Proof: Any Gaussian quantum channel can be written 
as £g(q) = ^E[Ub(g <E> te)U^], where £/& is a unitary 
operation corresponding to an Hamiltonian at most bi- 
linear in the field modes and where te is a Gaussian state 
[107| . Then, by using lemmas B2, B3 and B4 we obtain 

~6b[£g(q)] < 6[U b (g®T B )Ut] = 6 B [g]. □■ 
In turn, this lemma provides a necessary condition for 
a channel to be Gaussian: given a quantum channel 
£, and a generic quantum state g, if the inequality 
8b[£(q)] < 8b[q] is not fulfilled, the channel is nG. It 
is also worth to notice that the monotonicity is fulfilled 
only for a proper completely positive (CP) map. As we 
will see in Sec. I VIII even if we consider the conditional 
evolution corresponding to a Gaussian measurement op- 
erator, the nG of the output states may increase. Indeed, 
in this case, we do not consider a full CP-map, but only 
one Krauss operator corresponding to the chosen (Gaus- 
sian) measurement operator. 

Lemma B7 For a set of states {gk} having the same first 
and second moments, then nG is a convex functional, that 
is 

s B^2pkQk] < ^2pkS B [gi], 

k k 

with J2kPk = 1- 

Proof: The states gk, having the same first and second 
moments, have the same reference Gaussian state r which 
in turn is the reference Gaussian state of the convex com- 
bination g = ~^2kPkQk- Since conditional entropy S(p||r) 
is a jointly convex functional respect to both states, we 
have S B [J2 k Pk0k] = S(J2 k PkQk\\T) < EfcPfc s (efcll r ) = 

Notice that, in general, 8b is not convex, as it may easily 
proved upon considering the convex combination of two 
Gaussian states with different parameters. 

Lemma B8 At fixed average number of photons N — 
(aJa), the maximum value of nG measured by 8b for sin- 
gle mode states is achieved by pure superpositions of Fock 
states = Sfe a k\n+lk) where n > 0, Ik > ifc_i+3 or 
Ik = 0, and with the constraint N = (det <r[v(N)])^ — i — 

n + J2k 

Proof: Since 5b[q] — S(r) — S(g) we have to maximize 
S(t) and, at the same time, minimize S(g). For a 
single-mode system the most general Gaussian state can 
be written as g G = D(a)S(C)v{rk)&(£)Di (a), D(a) 
being the displacement operator, S(Q the squeezing 
operator, a,( € (D, and v(nt) a thermal state with n* 
average number of photons. Displacement and squeezing 
applied to thermal states increase the overall energy, 
while entropy is an increasing monotonous function 
of the number of thermal photons n t and is invariant 
under unitary operations. Thus, at fixed energy, S(r) is 



maximized for r = v(N). Therefore, the state with the 
maximum amount of nG must be a pure state (in order 
to have 5(g) = 0) with the same CM er = (N + ±)I of 
the thermal state v(N). One can easily check now that 
the state with this property is the one indicated in the 
Lemma. One can also observe that by choosing n = N 
and Ik = 0, we obtain that Fock states \N) are maximum 
nG states at fixed energy. □ 

As it will be clear from the examples presented in the next 
sections the two nG measures measures induce different 
ordering on the set of quantum states, that is we may find 
a pair of states g\ and gi such that <Sa[pi] > ^[^2] and 
^b[0i] < <5b[£>2], or viceversa. One may conje cture that, 
as it happens for entanglement measures |109j , we indeed 
do not have a unique nG measure and that different mea- 
sures correspond to different operational meanings. As 
also remarked in Sec. IVIIIi an operational meaning for 
8b may be found in terms of information-theoretic quan- 
tities, while an operational meaning for 8a, besides its 
connection with the distance in the phase-space, is still 
missing. The two measures are connected e ach o ther by 
means of the inequality S(q\\t) > D^ s [(?,r] [iToj . which, 
in turn, implies the inequality 

6b[q]>S A [q]h[q]. (39) 
For pure states fl3"9"]) reduces to Ssig] > SA[g]- 

C. A measure of nG based on the Wehrl entropy 

A different measure of nG has been proposed in p78j | . 
based on the difference between the Wehrl entropies of 
the reference Gaussian state and the quantum state in 
exam: 

6c[q]=H w (t)-H w (q) (40) 

where 

H w (g) = - [ d 2 aQ e (a)log{irQ e (a)} (41) 

is the Wehrl entropy, i.e. the differential entropy of the 
normalized Q-Husimi function 

Q e {a) = -(0\D( a y g D(a)\0) 

7T 

of the state g. The quantity 8q owns reasonable prop- 
erties in the phase-space, which are inherited from those 
of the the Q-function. However, it lacks an operational 
meaning and turns to be not invariant under Gaussian 
unitary operations. In order to illustrate this behaviour 
we have (numerically) evaluated the non-Gaussianity 8c 
for Fock number states \n) subjected to squeezing. In 
Fig. [3]we show Sc[S(r)\n)(n\S(ry] as a function of r for 
different values of n. As it is apparent from the plot, the 
nG is neither constant nor monotone with the squeezing 
parameter r. 
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A. Fock states and superpositions 



0.4 
0.3 
0.2 




0.2 0.6 1. 1.4 1.8 



FIG. 3. The Wehrl entropy-based non-Gaussianity 

8c[S(r)\n)(n\S(ry] for squeezed Fock states as a function of 
the squeezing parameter r. From bottom to top the nG for 
n = 1,2,3,4. 



We consider single mode Fock states 
sitions of Fock states of the form 



and superpo- 



\ipnk) 



1 

71 



|n) + \n + k) 



for n > and k > 2. The reference Gaussian states are 
thermal states T n k = u(n + |) with n + -| average pho- 
tons. NG can be analytically evaluated for both measures 
obtaining (for k = and k > 2) 



1 



1 



1 



2n + k 



k + 1 

SB[ipnk] = h(n+ — — ) , 



(42) 
(43) 



where the overlap O nk = {tp n k\v{n + t)|V>nfc) is given by 



D. Non-Gaussianity of a quantum operation 

Once one has at disposal a good measure for the non- 
Gaussianity of quantum state this may be exploited to 
define a measure for the non-Gaussian character of a 
quantum operation. Let us denote by Q the whole set 
of Gaussian states. A convenient definition for the non- 
Gaussianity of a map £ reads as follows 



5[£] 



max<5[£(g)] 

geG 



where £(g) denotes the quantum state obtained after the 
evolution imposed by the map. Indeed, we have 5[£ g ] = 
iff £ is a Gaussian map £ g , i.e. a map which trans- 
forms any input Gaussian state into a Gaussian state. 
Other properties follow from those of the nG measures 
for states. 

Despite the conceptual simplicity of the above defini- 
tion, the evaluation of S[£] is, in general, a challenging 
task using either the HS-based or the QRE-based mea- 
sure of nG. As a consequence, it has not been used so far 
for a systematic classification of maps. 



V. NON-GAUSSIANITIES OF SPECIFIC 
FAMILIES OF QUANTUM STATES 

This section is devoted to a sort of zoology of non- 
Gaussianity, i.e. we will consider different families of rele- 
vant quantum states and evaluate their non-Gaussianitics 
5a and 5b- In this way, we analyze in some details the re- 
lationships between the two measures and illustrate their 
basic features also in connection with the analytical prop- 
erties of their density operators and the intuition coming 
from their phase-space quasi-distributions. 



Onk = - 



( n+ !)r>+k 



As it is apparent from Fig. 2] both measures increases 
with both n and k and are monotone functions of each 
other for this families of states, i.e. if <5a[pi] > £>a[p2] 
then <5b[pi] > 5b[q2\- As stated in Lemma B7 Fock 
states have the maximum nG at fixed number of photons 
according to the measure 5b- Though it has not been 
proved yet, we observe the same result for the HS-based 
5a in all the examples considered up to now. 



3.5 
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2.5 
2.0 
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0.36 0.4 0.44 0.48 M 

FIG. 4. QRE based nG 5b as a function of HS distance nG 
5a for Fock states \n) with n = 1, . . . , 15 (red circles) and 
for superpositions \4> n k) with n = 1, ... 15 and k — 3 (blue 
squares), k — 4 (green diamonds), and k = 5 (black triangles). 
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B. Mixtures of Fock states 

We now investigate the monotonicity of the two mea- 
sures for other one-parameter families of quantum states. 
In this case we consider mixtures of Fock states of the 
form 



Qd 



oo 

E 

n=0 



q n (X) \n)(n\ 



(44) 



where < q n (X) < 1, J2 n 1™W = 1 an ^ n x = 
^2 n nc ln(X) is the average photon number of the state. 
The reference Gaussian state for any diagonal mixtures 
of Fock states is a thermal state v{n\) with the same av- 
erage photon number. The non-Gaussianity can be thus 
written as 



$a[qd] = 2 



1 



Sb [qd] = h(n x + 1/2) + ^2 In lo § <!■ 



(45) 



Si 
1.4 



B 




1. 
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FIG. 5. (Color online) Left panel: QRE based nG 5b as a 
function of HS distance nG 5a for mixtures of Fock states 
qd- Red circles refer to a Poissonian distribution, whereas 
empty black circles are for the diagonal states obitained as 
partial trace of TMC; green squares and black diamonds are 
for mixtures coming from partial traces of photon-added and 
photon-subtracted two-mode squeezed vacuum respectively; 
blue and purple triangles correspond to T-distributions (A) 
and gi 4 ' (A) respectively. Right panel: QRE based nG 5b as a 
function of HS distance nG 5 a for a sample of 10 4 truncated 
random mixtures qh of Fock states. The red cloud refers to 
H = 10, green for H = 100 and black for H = 1000. 



where r„ = {n\v{n\)\n) are the matrix elements of the 
(thermal) Gaussian reference state. We have numer- 
ically evaluated the non-Gaussianities for several one- 
parameter families including the diagonal states obtained 
as partial traces of TMS, PSSV and PASV states (see 
Eq. and the discussion following Lemma B4) as 

well as diagonal states with Poissonian profile or given 
by a r-distribution of the form oc n k exp{— n/X}. 

Results are shown in the left panel of Fig. [SJ the two 
measures are monotone each other for all the considered 
families and the behaviour of ^b[^_d] vs. 5a[qd] is almost 
independent on the kind of states. These results suggest 
that a general relation between the two measures may ex- 
ist for mixtures of Fock states. However, so far we have 
not been able to prove it analytically starting from the 
expressions of ^[^d] and £b[£>_d] in Eqs. (f4"5j) . 

We have also considered (truncated) random mixtures 
of the form 



H 

Qh = ^2 Pn 

71=0 



ri)(n\ 



where H is the truncation dimension, and have numeri- 
cally evaluated the non-Gaussianities using Eqs. (|4"5)) for 
a sample of 10 4 states. Results are reported in the right 
panel of Fig. [5] and show that despite the large number 
of involved parameters (up to H = 1000) the two mea- 
sures are almost monotone each other. As H increases 
the distribution of the two measures concentrates around 
the typical values. 



C. Schrodinger cat states 

Let us now consider the two-parameter family of quan- 
tum states given by the Schrodinger cat-like states, that 



is superpositions of coherent states |a) = D(a)\0) and 
I - a), 



cos co a) 



sin © — a) 



^1 + sin(20) exp{-2a 2 } 



(46) 



For (f) — i 7r /4, \ips) reduces to the so-called odd and 
even Schrodinger cat states. Using the fact that the refer- 
ence Gaussian state is a displaced squeezed thermal state 
t s = D(C)S(r)u(N)S^(r)D'f(C), where the real param- 
eters C, r and N are analytical function of <fi and a we 
have evaluated the nG measures 5a and 5b for different 
values of a and <j). Results are shown in Fig. [51 

As it is apparent from the plots, upon varying the value 
of the parameters the two measures exhibit similar qual- 
itative behaviour. In particular, for low values of the 
amplitude {e.g. a = 0.5), we observe an asymmetric be- 
haviour with respect to <6: nG is almost zero for positive 
(j), while for <6 < one achieve high values of nG. By 
increasing the value of the amplitude, say a = 5, both 
measures become even function of <j). This can be under- 
stood by looking at the Wigner functions of even and odd 
Schrodinger cat states. In fact for low amplitudes, the 
even cat (4> = 7r/4) Wigner function is similar to a Gaus- 
sian state, in particular to the vacuum state, while for 
4> = — 7r/4 it presents a non- Gaussian hole in the origin 
of the phase space; for higher values of a one can observe 
similar non-Gaussian beahviours both for the even and 
the odd cat state. 

Although the two nG measures capture the same qual- 
itative non-Gaussian behaviour, it is apparent from the 
parametric plot of Fig. [7] that they induce different or- 
dering on the set of states. In fact, upon varying the two 
parameters a and 4> one may find pair of states g\ , Q2 for 
which 5a [qi] > 5a[q2] and 5b [gi] < 5b[q2\- As discussed 
before, we do accept that the two measures may induce 
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FIG. 7. QRE based nG 8 b as a function of HS distance nG 8 a 
for Schrodinger cat states \4>s}- Solid blue lines refer to fixed 
amplitude of coherent states a, varying the angle —tt/2 < <f> < 
it/2; from bottom to top we have a — 0.5, 2.5. The dashed 
red lines are for fixed angles (j> and varying the amplitude < 
a < 2.5; from bottom to top </> — — 7r/3, 7t/6, 2n/5. The gray 
area denotes all the allowed values for the two nG measures 
for the considered range of values of the two parameters. 



FIG. 6. HS based nG 8a (left upper panel) and QRE based 
nG 8b (right upper right) for Schrodinger cat states \il>s) as a 
function of <f) and for different values of the amplitude a. Solid 
bue line: a = 0.5. Red dashed line: a — 5. The corresponding 
Wigner functions for the odd and the even cats are also shown 
in the lower panels: a = 0.5, (j> = +7r/4 (top left); a — 0.5, 
4> — — 7r/4 (top right); a — 5, cj> — +n/4 (bottom left); a = 5, 
4> — — 7r/4 (bottom right). 



different ordering on quantum states. Notice, however, 
that upon fixing one of the parameters and varying the 
other, one observes again a monotonous behaviour: the 
red dashed lines in Fig. [7] show 8b vs 8a for fixed values 
of <j) and varying the amplitude a. This appears to be 
a typical behaviour: for all the one-parameter families 
considered, the two measures are monotone with respect 
to each other and induce the same ordering of nG. 



VI. GAUSSIFICATION AND 
DE-GAUSSIFICATION PROCESSES 

In this section we will consider a single-mode Gaussi- 
fication, the loss mechanism due to the interaction with 
a bath of harmonic oscillators at zero temperature, and 
a de-Gaussification processes due either to phase diffu- 
sion or Kerr interaction. Other Gaussification and de- 
Gaussification protocols are considered in the next sec- 
tion where we discuss the case of entanglement distilla- 
tion. 



A. Loss mechanism 

The evolution of a single-mode quantum state inter- 
acting with a bath of harmonic oscillators at zero tem- 
perature is described by the following Lindblad Master 
equation 

g=^C[a]g (47) 

where g denotes time derivative, 7 is the damping factor 
and the superoperator C[0] acts as follows 

C[0]q = 20*00 -O^Oq- qO^O. (48) 

Upon writing rj = e~ 7 ' the solution of the Master equa- 
tion can be written as 

Q(r)) = J EV mQ Vl (49) 

m 

V ml 

where g is the initial state. In particular if the system 
is initially prepared in a Fock state g p = \p)(p\, after the 
evolution we obtain the mixed state 

p 

Q P (V) = E y ™^™ = Z>,pfo)IOv?l (50) 

m 1=0 

with 

^(i) = W(i-rV- (5i) 
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Since the state is diagonal in the Fock basis, the reference 
Gaussian state is a thermal state T p (rf) = v{prf) with 
average photon number prj. Non-Gaussianity Sa can be 
evaluated analytically 



Sa = 



1 



2{l-r)) 2m 2 Fi (-m,-m,l 



, < {l-i 1 ) 2m 2 F 1 [ -,,,.-,„. \ 

+ (1 + 2mr?)~ 1 - 



2(1 + (to - l)r])' 

(1 + TO77) m+1 



where 2F\(a,b,c,;x) denotes the Hypergeometric func- 
tion, and while 8 b can be evaluated numerically via the 
formula 



8b = PV log 



PV 



1 



PV 



log(l +pv) 



■^aiAv) !og(a;, P (?7))- 



1=0 



Because of Lemma B6 we know for sure that 8b is de- 
creasing with time r]t, while this property is not guaran- 
teed for 8 a- In Fig. [8]we plot both non-Gaussianities as a 
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FIG. 8. HS based nG 8 A (left panel) and QRE based nG 8 B 
(right panel) for a Fock state \p) under loss mechanism as a 
function of dimensionless time rjt and for different values of 
p. From bottom to top p — {2, 4, 6, 8}. 

function of dimensionless time rjt and for different values 
of p, observing that also 8 a is decreasing with time, and 
that both are monotonically increasing function of p, that 
is, at fixed time t the higher is the initial photon number 
p, the larger is the nG of the evolved state. Although 
they present a different shape, due also to their different 
scale (8 a for single- mode states is bounded by 1/2, while 
8b is in general unbounded), we observe a similar trend 
for both nGs which in particular approach zero for the 
same values of the parameters rj, t and p. 



B. Phase-diffusion evolution 

Let us consider single-mode systems evolving according 
to the following Master equation 



V C[a)a]g, 



(52) 



where the super-operator C[A] has been defined in Eq. 
(|4"5| and r is the noise factor. This Master equation 
describes the evolution of a quantum state subjected 
to a phase-diffusive noise. These non-Gaussian fluctu- 
ations are an important source of noise in optical com- 
munication links, and protocols able to purify squeezing 
or distilling entanglement have been recently proposed 
[it! l48j . Upon writing g in the Fock state basis, the ME 
leads to differential equations for the matrix elements 
Qnm = (n\g\m), where g nm = -|T (ra- m) 2 g nm whose 
solutions read as follows 



Qnm (^) 



3 -A 2 (n-m) 1 



£W(0). 



(53) 



In the last equation we defined A 2 = rt/2 whereas 
f?nm(0) denote the matrix elements of the initial state. 
^From Eq. (f53|) it is clear that the off-diagonal elements of 
the density matrix are progressively destroyed, whereas 
the diagonal ones are left unchanged and, in turn, energy 
is conserved. 

It is worth noticing that the same evolution as in (|53[) 
can be also obtained by the application of a random, 
zero-mean Gaussian-distributed phase-shift to the quan- 
tum state. Since the phase shift of an amount ip is de- 
scribed by the unitary operator U v = exp(— iip a'a), we 
can write the state degraded by the Gaussian phase noise 
as follows: 



QGn = dtp 



V4ttA 2 



U v g(0)Ul 



(54) 



nm JJR 



— iip(n—m) 



Vi^rA 2 



(0)\n)(m\ 



Qnm{0)\n){m\ . 



(55) 
(56) 



The parameter A is related to the width of the Gaussian 
distribution of the random phase-shift by the relation 
<y 2 rnd = 2A 2 : as one may expect, the more the Gaus- 
sian distribution is broad, the higher is the phase-noise 
affecting the quantum state. 

If we consider the mode initially prepared in a coherent 
state | a) with real amplitude the phase-diffused state is 
a non-Gaussian mixed state with density operator given 
by 

l2 ~ a n+m -A 2 (n-m) 2 

g A , a = e-^Y- %= \n)(m\, (57) 



n,m 



vnTro! 



in the Fock basis. 

We have evaluated numerically the nG <5s[pa,q] f° r dif- 
ferent values of the noise parameter A and of the average 
number of photons \a\ 2 . The results are shown in Fig. 
[5] In the left panel we report 8b as a function of A and 
for different values of a: large values of nG are achieved 
and the more intense is the initial coherent state, the 
more non-Gaussian is the output. For large values of the 
noise parameter the off-diagonal elements of the density 
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FIG. 9. (Left panel): QRE based nG <5g for coherent states 
undergoing phase-diffusion as a function of the noise param- 
eter A and for different values of the average number of pho- 
tons. From bottom to top: \a\ 2 — {1,2,3,4,5}. (Right 
panel): QRE based nG Sb for coherent states undergoing 
phase-diffusion as a function of the average number of pho- 
tons \a\ 2 for different values noise parameter: dashed lines, 
from bottom to top: A = {0.25, 0.5, oo}. The black solid 
line represent the maximum value of nG at fixed number of 
photons. 



matrix are completely destroyed and the nG approaches 
its aysmptotic value, corresponding to the nG of a diago- 
nal mixture of Fock states with a Poissonian distribution. 
The asymptotic value may be evaluated using Eq. (j45|) 
and choosing q n as a Poissonian distribution with mean 
value A = \a\ 2 . Before approaching the asymptotic value 
$b[qa,o\ is not monotone: a not much pronounced maxi- 
mum of Sb may be seen for intermediate values of A. In 
the right panel of Fig. [HI we plot the nG for different val- 
ues of A as a function of the number of photons |a| 2 . As 
we noticed above the maximum value of Sb [pa,q] at fixed 
\a\ 2 is not the asymptotic value. However, the difference 
is very small and thus, upon observing the behaviour for 
different values of A, we conclude that the maximum nG 
cannot be achieved by this family of quantum states. 



C. Kerr interaction 

One of the simplest unitary non-Gaussian evolution is 
provided by the so-called self-Kerr effect taking place in 
third-order nonlinear media. The interaction Hamil- 
tonian is given by 



{a)af 



(58) 



and the evolution operator by Uk. er r — exp{— (era) 2 } 
where 7 = Tt is a dimensionless coupling constant. Kerr 
interaction has been suggested to realize quantum non- 
demolition measurements, to enhance quantum estima- 
tion performances in qu antum op tics and to generate 



E iuan tum superpositions [111 
1131 ] and entanglement [5£ 



|112j as well as squeezing 
i9]. A known example of 



Kerr medium is provided by optical fibers where, how- 
ever, nonlincarities arc small and accompanied by other 
unwanted effects. Recently, larger Kerr nonlinearities 
have benn proposed in many different physical systems 
[l!4| and have been observed with electro-magnetically 



induced transparency [115J , with Bos e-Einstein conden- 
sates jll6l | and with cold atoms |l!7j . These results re- 
newed the interest for the quantum effects of Kerr inter- 
action, which are always accompanied by the generation 
of non-Gaussianity. In the following, we consider the nG 
features of an initial coherent state \a) undergoing Kerr 
interaction 



Ukerr\ot) 



\a\ 2 /2 ST 



|n). (59) 



Since the evolution is unitary, the output state is still 
pure and the evaluation of the nG Sb is straightforward 
and involves only the computation of the covariance ma- 
trix of the evolved state. Non-Gaussianity <5b[|<Xv)] is 



n 



as a function of the average num- 

v2| 



plotted in Fig. 

ber of photons n = |a 2 | and for different values of the 
coupling constant 7. The maximum nG achievable for 
given number of photons is also reported for comparison. 
As expected, non-Gaussianity is an increasing function 
of the initial energy n and, for the range of values here 
considered, of the the coupling constant 7. 




n 



FIG. 10. QRE based nG Sb of coherent states undergo- 
ing Kerr interaction as a function of the average number of 
photons and for different values of the coupling constant 7. 
Dashed lines from bottom to top 7 = {10~ 6 , 10" 4 , 10" 2 }. The 
black solid line is the maximum nG at fixed number of pho- 
tons. 

For 7 « 10 -2 the maximum nG achievable at fixed 
energy is quite rapidly achieved, while for more realistic 
values of the Kerr coupling nG is obtained only for large 
values of the average number of photons. In the exper- 
iments proposed to obtain entanglement via Kerr inter- 
action [58l l59j. pulses with an average number up to 10 s 
photons are needed to compensate the small nonlineari- 
ties of standard glass fibers. Therefore, in these regimes, 
the generation of entanglement is always accompanied by 
a large degree of nG. 
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VII. NG AND DISTILLATION OF 
ENTANGLEMENT 



Long distance quantum communication requires that 
the communicating parties share highly entangled states 
over long distance. One has therefore to deal with the 
daunting task of distributing highly entangled states over 
long distances, overcoming losses and decoherence due to 
the unavoidable coupling of the system with the environ- 
ment. For discrete-variable states (qudit) entanglement 
distillation protocols, where a subset of states with an 
high degree of entanglement are extracted from an en- 
emble of less entangled states, have been proposed and 
experimentally demonstrated. As regards continuous- 
variable entanglement, it has been proved that entangle- 
ment distillation cannot be performed within the Gaus- 
sian world, i.e. by starting from Gaussian entangled 
states and by linear optical components, homodyne de- 
tection and classical communication. Non-Gaussianity 
is a necessary ingredient in an entanglement distillation 
protocol. In particular, two different main approaches 
can be adopted: in the first approach, e.g. the protocol of 
HI], one starts with an entangled non-Gaussian state, 
and then use Gaussian operations based on linear opti- 
cal elements, homodyne detection or vacuum projective 
measurements. In the second approach one starts with an 
entangled Gaussian beam and try to increase its entangle- 
ment by using non-Gaussian operations such as photon 
number conditional measurements. In [25— 27| it has been 
proved that two-mode squeezed states with photons sub- 
tracted on the two modes can be used to obtain better 
teleportation fidelities. A first full scheme of entangle- 
ment distillation based on this idea has been presented 
in 49] , where increase of entanglement by means of lo- 
cal photons subtraction from two-mode Gaussian states 
has been observed. In the following we will review the 
protocols of [13] and pointing out the role played 
by non-Gaussianity and its amount in the success of the 
protocols. 

The protocol proposed in [47j . from now one the B- 
protocol, makes use of beam splitters and on/off detec- 
tors, i.e. detectors only able to distinguish the presence 
or the absence of photons. The input state of the protocol 
is the state g <g> g i.e. two replicas of a two-mode non- 
Gaussian. The two copies are mixed in a balanced beam 
splitter and then two of the output modes are directed 
into on/off photon detectors. The state is kept if both 
the local detectors are registering the outcome "zero", 
i.e., at least in ideal conditions, the presence of the vac- 
cum state. In [13] it has been proved that the B-protocol 
drives the initial state towards a zero-displacement Gaus- 
sian state, which in turn are the only fixed points of the 
map. Moreover, under some assumptions on the initial 
state, the protocol acts as a purifying protocol as well, 
distilling the initial non-Gaussian entangled mixed state 
into a pure Gaussian entangled state. Notice that despite 
the fact that the output state of each step is obtained 
via (Gaussian) passive linear operations and Gaussian 



measurements (the projection on the vacuum state), the 
map cannot be descr ibed by a proper completely positive 
Gaussian operation [l 18j | . 
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FIG. 11. Non-Gaussianity in the B-protocol. The plot shows 
the QRE based non-Gaussianity 8b of the output state after 
s steps of the B-protocol, when the initial state g a is given 
in Eq. ()6Up . as a function of the parameter A. Step: s = - 
black solid line; s = 5 - black dashed line; s = 10 - gray solid 
line; s — 20 - gray dashed line. 



As an illustrative example let us consider the B- 
protocol applied to the non-Gaussian pure state with 
the following non-zero matrix elements in the Fock basis 
g{a, b, c, d) = (a, b\g\c, d) 



Qa(0, 0,0,0) 



1 



1 + A 2 



g a (0, 0,1,1) = g (l,l, 0,0) = 



1 + A 2 



0a(l, 1,1,1) 



A 2 



1 + A 2 



(60) 



In Fig. [IT] we show the non-Gaussianity quantified by 5 b 
as a function of the parameter A. for states obtained after 
different number of steps of the protocol. As a matter of 
fact, at fixed A, nG is not always monotonically decreas- 
ing under the iteration of the protocol, and for A 1 the 
value of nG may increase, even achieving very high val- 
ues. On the other hand, the overall effectiveness of the 
protocol is confirmed by our measure, since the range of 
values of A for which 5 b ~ increases at each step of the 
protocol. Overall, the use of our nG measure may help 
to tailor the distillation protocol depending on the initial 
conditions. 

Let us now consider as initial states the pure state g a 
of Eq. I|60p and the mixed state gb with non-zero matrix 



15 



elements 



,(/) 



56(0,0,0,0) = 



1 



1 + A 2 

g 6 (0, 0,1,1) = 56(1, 1,0,0) 

A 2 

56(1,1,1,1) 



2(1 + A 2 ) 



1 + A 2 



(61) 



Both states converge towards pure Gaussian entangled 
states, and it has been shown [47j that entanglement in- 
creases at each step of the protocol. Here we will inves- 
tigate how much the gained entanglement is related to 
the non-Gaussianity of the initial state. For this pur- 
pose, since both, entanglement and non-Gaussianity are 
increasing quantities with the number of photons, we will 
consider a renormalized version of non-Gaussianity 5r[q] 
and the relative entanglement gain at each step A W . The 
maximum amount of non-Gaussianity for a two-mode 
state with N = Tr[g(cva + b'b)] photons is 



2[(1 + N/2) log(l + N/2) - (N/2) log(JV/2)], 



and the renormalized non-Gaussianity is defined as 

5b[q] 



S ( m\n) 



(62) 



We define the relative entanglement gain at the step i as 

(63) 



where 5M is the output state at the i-th step of the pro- 
tocol. The degree of entanglement is quantified in terms 
of the logarithmic negativity, i.e. 



E N (g) = log 2 |j5 r ||i 



(64) 



where || ||i de notes the trace- norm, and g T is the partial 
transpose of g [l!9j . In Fig. [T2]we plot for both states g a 
and 5b, the increase of entanglement AW as a function 
of the renormalized non-Gaussianity <5r[5]. We observe 
that the more the initial state is non-Gaussian, the larger 
is the entanglement increase at each step of the protocol. 
Therefore, we may say that at least for this particular 
protocol non-Gaussianity plays a relevant role, and it is 
quantitatively responsible for the good performances of 
the protocol. 

Let us now consider the distillation protocol proposed 
in 49], from now on the T-protocol. In this protocol, an 
entangled Gaussian state is obtained by mixing a single- 
mode squeezed state \tp) = S(r)\0) with the vacuum in 
a balanced beam-splitter. Then, both the modes are 
splitted to perform a photon-subtraction and the state 
is keeped if one or both detectors clicks (see [49| for de- 
tails). The two possible output states of the distillation 
procedure can be written in a compact form as 

\^lt)AB =Ma n A A a n B B B AB (n/A)S A (r)\0) A \0) B (65) 
= NB AB (tt/4) a\ A +nB S A (r) |0> A \0) B (66) 
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FIG. 12. The relative increase of entanglement at the 
i-th step of the B-protocol as a function of the renormalized 
nG of the initial state (left panel: qa\ right panel: qb)- From 
bottom to top: i = {1, 2, 5, 00} 



where if (n A ,n B ) = (1,0) = (0,1) we have a single- 
photon subtracted state , while for (n A ,n B ) = 
(1,1) we obtain the two-photon subtracted state \ip ut>- 
The non-Gaussianities of the two states can be easily 
evaluated by using the properties of the measure S B , and 
in particular by exploiting the invariance under Gaus- 
sian unitary operation such as beam-splitter evolution 
and squeezing. After some calculations, it turns out that 



5 



= 8 B \\l) 



6BU { ol)}=5 B [Af'UA0)+V2v\2))} 

where fi = cosh(r), v — sinh(r) and M' is a normalization 

factor. We can immediately observe that for |V>o«t) the 
nG does not depend on the squeezing parameter r. In 
Fig. Q2] we plot as a function of r respectively in the 
left panel the entanglement and in the right panel the 
non-Gaussianities of the two possible ouput states. We 
observe that for small values of r the entanglement and 
the nG of the output states have a similar beahviour. 
However this similarity is lost when we consider higher 
values of the squeezing, i.e. we can say nothing in general 
about the relationship between the distilled entanglement 
and the non-Gaussianity for this protocol. 
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FIG. 13. Entanglement (left panel) and QRE based non- 
Gaussianity (right panel) of the distilled states in the T- 
protocol as a function of the squeezing parameter r . Dotted 
red: one-photon subtracted state 
photon subtracted state \ip^t)- 



V^Jt). Solid blue: two- 



In summary, the role of the amount of non-Gaussianity 
in continuous variable entanglement distillation protocols 
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is still not fully understood and worth to study for pro- 
tocols employing non-Gaussian entangled states as ini- 
tial resource. In particular, an interesting open question 
that arises in this framework is whether there is a max- 
imum amount of distillable entanglement at fixed non- 
Gaussianity of the resource state. 

VIII. NG AND QUANTUM COMMUNICATION 

As anticipated in the introduction, Gaussian states 
play a relevant role in quantum communication, in par- 
ticular for being extremal, at fixed covariance matrix for 
several relevant quantities [|J. For example, channel ca- 
pacities are maximized by Gaussian states, whereas most 
of the entanglement measures are minimised by Gaus- 
sian states. In this section we will prove the extremality 
conditions for some of these quantities by using our non- 
Gaussianity measure 8b, and studying the role played by 
non-Gaussianity itself in these inequalities 

As a first example, let us consider a very important and 
well known quantity in quantum information, i.e. the x~ 
Holevo quantity. Let us suppose to transmit information 
by encoding the symbols a^, chosen with probabilities pi 
from an alphabet A, in a set of quantum states Qi. At 
each use of channel, the state preparation is described by 
the overall state g = J2iPi°i> and the Holevo quantity, 
which quantifies the amount of accessible information, is 
given by the formula xio) — — ^2sPiS(Qi)- If we fix 
the covariance matrix of the overall state g and consider 
pure encoding states Qi, we obtain (we will denote 5 b = 5 
in the rest of the paper), 

X (g) = S(r) - S[g] (67) 

The maximum value of x is attained by considering the 
Gaussian state and the nG S[g] exactly quantifies how 
much information is lost by considering a non-Gaussian 
overall state. 

Let us now consider bipartite states and rephrase the 
extremality condition for two quantities: the quantum 
mutual information and the quantum conditional en- 
tropy. The quantum mutual information, given a bipar- 
tite stat e pa b is defined as I(A : B) = S(qa) + $(qb) — 
5(qab) [l 201 ] . It quantifies the amount of correlations 
(classical and quantum) in a bipartite state a nd, f or pure 
states corresponds to the entanglement. In |l2l| it has 
been shown that at fixed covariance matrix I (A : B) is 
minimised by Gaussian states. By using our measure 5 
the proof is simply based on the lemma B5, that is 

I{A:B) = S{qa) + S{q b )-S{qab) 
= S(t a ) + S{t b ) - S(t ab ) 

+ {8[g A B} - 8[g A ] ~ 8[g B }) 
> Ig{A : B) (68) 

where Ig(A : B) = S(ta)+S(tb) — S(tab) is the quantum 
mutual information obtained by considering the reference 



Gaussian states. The two mutual informations I (A : B) 
and Ig(A : B) differs exactly by the quantity 

A 2 = 5[qab] - 8[g A ® qb] , 

that is, the amount of correlations that are lost upon 
considering the Gaussian counterpart of gAB is equal to 
the amount of non-Gaussianity that is lost by considering 
the tensor product of the partial states gA <X> £»b instead 
of the (correlated) bipartite state gAB ■ 

Another quantity that is known to be maximised, at 
fixed covariance matr ix, by Gaussian states is the con- 
ditional entropy [107| . Conditional entropy is defined, 
given a bipartite state gAB as S(A|-B) = S(qab) — $(Qb) 
|l20j | . The proof of the extremality can be easily obtained 
by means of our QRE based measure in the following way: 

S(A\B) = S(qab)-S(q b ) 

= S(t A b) - S(t a ) - {8[qab\ - 8[qb\) 

< S G (A\B) (69) 

where we used the monotonicity of the measure under 
partial trace, and where we defined the Gaussian condi- 
tional entropy Sg{A\B) = S(tab) — S(t~b) as the one eval- 
uated by considering the reference Gaussian states. The 
two conditional entropies differ exactly by the quantity 
Ai = 5[^ab] ^8[gs]\ thus, the more the non-Gaussianity 
of the overall state is robust under discarding a subsys- 
tem, the more the difference between the two quantities 
Sg(A\B) and S(A\B) is near to zero. This extremality 
condition is important for several reasons: the negative 
of the conditional entropy is a lower bound for the dis- 
tillable entanglement, and thus, one can evaluate, given 
a generic bipartite state qab, a simple lower bound on 
the distillable entanglement based only on first and sec- 
ond moments of the state. Moreover, according to the 
op erati onal meaning of quantum relative entropy given 
in [l22t | in the framework of quantum state merging, we 
observe that at fixed covariance matrix, it is always more 
convenient to use non-Gaussian state. Indeed, for posi- 
tive values, the quantum conditional entropy quantifies 
how much quantum information Alice needs to send to 
Bob so that he gains the full knowledge of the bipartite 
state gAB given his previous knowledge about the partial 
state gB, while for negative values it turns out that Al- 
ice needs to send only classical information and moreover 
the two users gain entanglement to perform, as example, 
teleportation. As a consequence at fixed covariance ma- 
trix, by using a non-Gaussian state you have to send less 
information or, for negative values, you gain more en- 
tanglement; the gain is exactly given by the degradation 
of the non-Gaussianity under partial trace, that is how 
much non-Gaussianity is lost from the initial state gAB 
tracing out the first Hilbert space. 

A different interpreta tion of the conditional entropy 
has been also given in jl23| . in the context of the so- 
called private quantum decoupling. It has been shown 
that, in the limit of infinitely many copies of the initial 
state gAB , the ineliminable correlations between the two 
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parties are quantified by the negative of the conditional 
entropy. Again we have that, at fixed covariance matrix, 
non-Gaussian state have more ineliminable quantum cor- 
relations than the corresponding Gaussian one, and that 
the more nG is lost under partial trace operation, the 
more these quantum correlations are present. 

The above results lead to speculate about non- 
Gaussian correlations which, being encoded in a larger 
set of degrees of freedom, are, at fixed covariance matrix, 
higher than the Gaussian ones. In particular, we have 
observed this non-trivial connection between the robust- 
ness of nG under partial trace and decoupling operations, 
as well as in the quantum correlations present in non- 
Gaussian bipartite states. Recently, arguments have been 
provided supporting the conjecture that at fixed energy 
Gaussian entanglement is the most robu st against noise 
in a Markovian Gaussian channel |l24j . On the other 
hand, the same analysis have also shown that robustness 
of non-Gaussian states is comparable with that of Gaus- 
sian states for sufficiently high energy of the states. This 
implies that in these regimes non-Gaussian resources can 
be exploited to improve quantum communication proto- 
cols approximately over the same distances. 



IX. NG AND PARAMETER ESTIMATION 

In an estimation problem one tries to infer the value 
the of a parameter A by measuring a different quantity 
X, which is somehow related to A. This often hap- 
pens in quantum mechanics and quantum information 
where ma ny quantities of intereset, e.g. entanglement 
[l25l Il26j |. does not correspond to a proper observable 
and should be estimated fro m the measurement of one 
or more observable quantities jl27l |. Given a set {gx} of 
quantum states parametrized by the value of the quan- 
tity of interest, an estimator A for A is a real function 
of the outcomes of the measureme nts perfor med on Q\. 
The quantum Cramer- Rao theorem jl28l . Il 29f | establishes 
a lower bound for the variance Var(A) of any unbiased es- 
timator, i.e. for the estimation precision, 



Var(A) > 



1 



MH{\) 



(70) 



in terms of the number of measurements M and the so- 
called quantum Fisher information, which captures the 
statistical distinguishability of the states within the set 
and itself is proportional to the Bures distance between 
states corresponding to infinitesimally close values of the 
parameter, i.e. 



H(X) = 4:d%(Qx + d\,Q\) 



\(lp m \dxQx\tp n )\' 

nm 



Qn + Q-n 



(71) 



where we have used the eigenbasis g\ = J2 n Qn\4'n){'4'n\- 
In an estimation problem where the variation of a pa- 
rameter affects the Gaussian character of the involved 



states one may expect the amount of non-Gaussianity to 
play a role in determining the estimation precision. This 
indeed the case: the non-Gaussianity 5b provides an up- 
per bound to the quantum Fisher information at fixed 
covariance matrix. This is more precisely expressed by 
the following 

Theorem 1 // t\ is a Gaussian state and an infinites- 
imal variation of the value of A drives it into a state 
Q\+d\ with the same covariance matrix, then the non- 
Gaussianity 5s[gx+dx] provides an upper bound to the 
quantum Fisher information. 

Proof: If gx+dx and tx have the same CM then the nG 
of gx+dx, S B [gx+dx] = S{gx+dx\\Tx) equals t he so - called 
Kubo-Mori-Bogolubov information H(X) \)M, EH, 
which itself provides an upper bound for the quantum 
Fisher information H(X) < H(X) [l32T |. thus proving the 
theorem. □ 

The theorem says that the more non-Gaussian is the 
perturbed state, the more may be distinguishable from 
the original one, thus allowing a more precise estimation. 
One may wonder that when gx+dx is itself a Gaussian 
state the theorem requires H(X) = 0, i.e. no reliable 
estimation is possible. Indeed, this should be the case, 
since Gaussian states are uniquely determined by the first 
two moments and thus the requirement of having the 
same covariance matrix implies that Tx+dx and tx are 
actually the same quantum state. 

For situations where the CM is changed by the pertur- 
bation we have no general results. On the other hand, 
it has been already shown that non-Gaussian states im- 
prove quantum estimation of displacement and squeezing 
parameters [3l[ and of the loss parameter [32], compared 
to optimal Gaussian probes. 



X. EXPERIMENTALLY FRIENDLY LOWER 
BOUNDS TO QRE NG 

A drawback of the nG measure 5b is that its evalua- 
tion requires the knowledge of the full density matrix g. 
For this reason, it is often hard to compute when one has 
only partial informations coming from some, maybe inef- 
ficient, measurements. In literature differente approaches 
have been proposed to estimate squeezing [1331 ] and en- 
tanglement [80( of Gaussian and non-Gaussian states 
when only certain measurement are available in the lab. 
In the following we will derive some lower bounds to the 
QRE-based nG measure for some class of states and by 
considering the possibility to perform on the states only 
certain efficient or inefficient measurements. 



A. Diagonal state and inefficient photodetection 

Let us consider a generic single-mode state diagonal in 
the Fock basis g = ^2 n Pn\n)(n\. Its non-Gaussianity can 
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be evaluated as 

S[q] = S(v N ) - S(g) = SK) - H(p n ) (72) 

where vn is a thermal state with the same average photon 
number N = n Pm arL d from now on 

H(p n ) = ~ ^Pn l0gp„ 

n 

denotes the Shannon entropy corresponding to the dis- 
tribution {p n }- 

Let us consider now an inefficient photodetector de- 
scribed by POVM operators 

oo 

n m = a m ,s (v) I s )( s I 

with a mtS {f]) defined in Eq. (|5ip . and where rj is the 
efficiency of the detector. Using this kind of detection 
one can reconstruct the probability distribution q m — 
Tr[f?II m ]. We want to show that the quantity 



€ A [q] = S(v M ) - H(q m ) 



(73) 



with M = mq m = rjN, is a lower bound on the ac- 
tual non-Gaussianity S[g\. To this aim we remind that 
an inefficient photodetection can be described by mix- 
ing the quantum state g with a vacuum state at a beam 
splitter with transmissivity r\ followed by a perfect pho- 
todetection with projective operators P m = \m){m\. The 
corresponding probability distribution is therefore 

q m = Tr 12 [U BS (ri)Q® |0)<0|t4 s (»j)|m)H <8 1] 
= Tr 1 [£(g)\m)(m\) 

where £{q) = Ty 2 [U B s(v)q ® |0)<0|?7^ g (r?)] is the loss 
channel applied on the quantum state g. Since g is di- 
agonal in the Fock basis we can easily show that £ (g) is 
still diagonal, 

n 

£(0) = 5>„£(|n)(n|) = £EP»«i.«(*7)|i><*l (74) 

n n 1=0 

where we used that £(\n)(n\) = Ya=o a l,n(vW){l\- By 
performing an efficient photodetection we can experimen- 
tally obtain the probability distribution 



B. State with a thermal reference Gaussian state 
and ideal photodetection 

Let us consider a quantum state g = J2 n mP n , m \ n )( m \ 
such that Tr[ga] = Tr[ga 2 ] = 0. The corresponding refer- 
ence Gaussian state is a thermal state i/jv with the same 
average number of photons N — Tr [ga^a]. Let us con- 
sider now the quantum state g\d obtained considering 
only the photon number distribution p n>n and removing 
all the off-diagonal elements that is, g\d — ^2 n Pn,n\n){n\. 
The reference Gaussian state of g\d is again the thermal 
state vn- Let us consider now the state 



A/a(p) = £ 



-A'(n-m)' 
e Pn,' 



\n)(m\ 



(76) 



which physically corresponds to a phase-diffusion applied 
to the initial state. As a matter of fact, the same kind 
of evolved state may be obtained by the application of a 
random zero-mean Gaussian distributed phase-shift to g, 
that is: 



^A(g) = 



-<j> 2 /(4A 2 



U*gUl, (77) 



where U<f, = exp{— ia^acf)}. Upon using the invariance 
under unitary operators and the concavity of the von 
Neumann entropy one may show that 



sCA/kte)) > 



-<p 2 /(4A 2 



■S(U 4> gUl)=S(g). (78) 



ir V4ttA 2 
^From Eq. (|T6[) we have that g\d = limA->co J^a(q) an d 



because of Eq. (f?8j) 

H(p„, n ) = S(g\ d ) > S(f>) 



(79) 



It is straightforward to see that the non-Gaussianity of 
g is lower bounded by the non-Gaussianity evaluated by 
considering only the photon-number distribution of the 
state, i.e. 

e B [g] = S(v N ) - H( Pn . n ) < S(u N ) - S{g) = S[g] (80) 



State with a thermal reference Gaussian state 
and inefficient photodetection 



q m = Tr[£ (g)\m)(m\] = ^ PnUm^in) 

The quantum state £(g) in fully described by q m and, 
by observing Eq. fl72"j) and (|73|) we can easily see that 
e A[g] = S[£(g)]. By simply using the fact the non- 
Gaussianity measure 5[g] is non-increasing under Gaus- 
sian maps, such as £, we finally obtain 



e A [g] = S[£(g)} < 5[g}. 



(75) 



This inequality tells us that performing an inefficient pho- 
todetection on a given quantum state g diagonal in the 
Fock basis and if we are able to reconstruct the proba- 
bility distribution q m , we can use Eq. (|73p to obtain a 
lower bound on the actual non-Gaussianity 5[g\. 



Let us take a quantum state g as in the previous Sec- 
tion, but considering the case of an inefficient photode- 
tection . ^From the measurement, we obtain the distri- 
bution q m — Tr[gTl m ], and we can thus define a quantum 
state 



q m \m)(m\. (81) 



having M = Yl m m 1 m avera g e photons. As in Sec. IX Al 
we have Qm = £{g\d), with g\d defined as before and £ 
denoting the loss channel. Then, by using the mono- 
tonicity of non-Gaussianity under Gaussian maps, and 
the previous results, we obtain 



ec[g] = S(v M ) - H(q„ 



hi] < S[g\ d ] < S[g] (82) 
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D. Generic state with known covariance matrix 

and ideal photodetection 

Let us consider a generic single-mode state g — 
mP n , m \ n )( m \ w ith a reference Gaussian state r. Its 
non-Gaussianity is evaluated as 

5[g] = S(r) - S(g) (83) 

Let us consider the quantum state g\d — ^2 n Pn,n\n){n\ 
as in the previous case and suppose that we able to eval- 
uate its covariance matrix of g (and thus the entropy of 
the reference Gaussian state r). We proved before that 
H(p n ,n) = S(g|d) > S(g) and thus we have the com- 
putable lower bound on the non-Gaussianity 

e D [Q\=S(T)-H(p n , n )<S[g] (84) 

E. Generic state with inefficient photodetection 

Let us consider a generic single-mode state g — 
J2 n m P n , m \ n }( m \ with a reference Gaussian state r. Be- 
cause of the monotonicity of the measure under Gaussian 
maps, we have 

5[g] > 6[£(q)} = S(r„) - S(£{q)) (85) 

where t v = £(r). Again, by using the inequality derived 
in Eq. (|79]) we obtain 

e E [g]=S(T v )-H(q m )<S[g] (86) 

where q m — Tr[glT. m ]. This general lower bound can be 
useful when the covariance matrix of the state can be 
easily derived from the photon number statistics of g (e.g. 
for phase-averaged coherent states). 

XI. CONCLUSIONS 

Non-Gaussianity is a resource for quantum information 
processing and thus we urge a measure able to quantify 
the non-Gaussian character of a quantum state. In this 
paper we have addressed non-Gaussianity of states and 
operations in continuous-variable systems and we have 
illustrated in details two measures of nG proposed in 
[76l [77} , along with their properties and the relationships 
between them. We used them to assess some Gaussifica- 
tion and de-Gaussification processes, and in particular we 
studied the role of the amount of nG in two entanglement 
distillation protocols proposed in literature. The role of 
non-Gaussianity appears to depend on the protocol it- 
self, and at least in one of the two protocols, the amount 
of gained entanglement at each step of the protocol is 
monotonous with the nG of the initial low-entangled 
state. We have also reconsidered the extremality of Gaus- 
sian states in terms of our measure based on the QRE, 
for some relevant quantities in quantum communication, 



as conditional entropy, mutual information and Holevo 
bound. In particular, we found that in the bipartite set- 
ting there is a, probably not entirely understood, connec- 
tion between correlations and nG: at fixed covariance ma- 
trix non-Gaussian states have more correlations and this 
excess of correlations is related on the amount of non-G 
that the quantum state loses under partial trace opera- 
tion or under decoupling. These results, together with 
recent ones on the robustness of non-Gaussian entangle- 
ment in noisy Markovian channels implies that there are 
regimes where non-Gaussian resources can be exploited 
to improve quantum communication protocols. We have 
also seen that QRE non-Gaussianity is a bound for the 
quantum Fisher information at fixed covariance matrix 
and thus the nG features of quantum states may be also 
used to improve parameter estimation with continouos 
variables. Finally, since the evaluation of the QRE nG 
measure requires the knowledge of the full density matrix, 
we derive some experimentally friendly lower bounds to 
nG for some class of states and by considering the pos- 
sibility to perform on the states only certain efficient or 
inefficient measurements. 

Our analysis of the properties and the applications of 
the two measures of non-Gaussianity has shown that they 
provide a suitable quantification of nG for the purposes 
of quantum information. In particular, the QRE based 
nG 5b have several operational characterizations and it 
is not too difficult to be evaluated. Since all states g 
with the same first two moments at fixed purity have the 
same amount of nG 5b[q] an interesting and remarkable 
picture emerges: for many purposes the effects of nG may 
be described by a single global parameter rather than 
ascribed to a specific higher moment. This is perhaps 
our main conclusion. Overall, our results suggest that in 
terms of resources for quantum information, the amount 
of non-Gaussianity of a quantum state can be evaluated 
using 5b, with 5a serving as a fine-tuning tool for specific 
purposes. 

There are several open problems requiring further in- 
vestigations about non-Gaussianity of quantum states. 
Among them we mention the following ones, which also 
provides a summary of the unanswered quastions posed 
in our paper: i) Is there any general relation between 
the two measures of non-Gaussianity, at least for specific 
class of states? ii) Is there a maximum amount of distill- 
able entanglement at fixed non-Gaussianity ? iii) Which 
is the role of non-Gaussianity in parameter estimation 
involving a change in the covariance matrix? 

In conclusion, non-Gaussianity is a resource that can 
be quantified. Our results pave the way for further de- 
velopment and suggest that a deeper understanding of 
the geometrical and analytical structures underlying the 
non-Gaussian features of states and operations could be 
a powerful tool for the effective implementation of quan- 
tum information processing with continuous variables. 
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